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1 Introduction 

The purpose of this article is to study the problem of constructing A^-dimensional, 
N > 2 , vectorial Sturm-Liouville differential equations subject to certain bound- 
ary conditions, which are isospectral to a given one. 

Those A^-dimensional vectorial Sturm-Liouville eigenvalue problems which are 
considered in this paper, are of the following form: 

(1) -(l)"{x) + P{x)^{x)=X^{x), B(t)'{0)+A(t){0) = B(l)'{Tr) + A(t){Tr) = O, 

where < x < vr , (^{x) is an R -valued function, P{x) is a continuous N x N 
symmetric matrix-valued function, A,B,A,B are N x N matrices which satisfy 
the following conditions: 

(2) BA* = AB\BA* = AB* , rank[^, B] = rank[^, B] = N, 

where A* is the transpose matrix of A, [A, B] denotes the N x 2N matrix whose 
first N X N block is A, and the second N x N block is B. We shall use the tuple 
{P, A, B,A, B) to denote the eigenvalue problem (|l|). Note that the conditions in 
(§) ensure the problem (|l|) a selfadjoint eigenvalue problem, and its eigenvalues 
can be determined by the variational principle. Counting multiplicities of the 
eigenvalues, we arrange the eigenvalues of (|l]) in an ascending sequence 

(3) /Uo < Ail < ^2 < • • • ■ 

This sequence shall be denote by S(P, A, B,A,B), and is called the sequence of 
eigenvalues of (|l]). Note that the multiplicity of each eigenvalue of (||) is at most 
A^. For convenience, we shall use a{P, A, B, A, B) to denote the set of eigenvalues 
of (H), and arrange its elements in ascending order as 

Ao < Ai < A2 < • • • , 

and use m^ denote the multiplicity of A^ in the sequence ^. Given two N- 
dimensional selfadjoint vectorial Sturm-Liouville eigenvalue problem (P, A, B, A, 
B) and {P,A,B,A,B) over < x < vr. li J:{P, A, B , A, B) = ^{P,A,B,A,B), 
we call these two eigenvalue problems isospectral problenas, or {P,A,B,A,B) 
is isospectral to {P,A,B,A,B). For scalar Sturm-Liouville equations, i.e., (||) 
with A^ = 1, isospectral problems have been studied by many mathematicians, 
notably, G. Borg ||l|, L M. Gel'fand , B. M. Levitan and their associates, and the 
structure of the set of isospectral problems for scalar Sturm-Liouville equations 
is well-presented in the book of J. Poschel and E. Trubowitz [g], and in the 
works of E. Trubowitz. But for vectorial Strum-Liouville equations, i.e., (|lj) with 
iV > 2 , methods for constructing isospectral problems and the structure of 



isospectral problems are not well-understood. Motivated by a recent work of 
Jodeit and Levitan y], in this paper we present a method for constructing an 
A'^-dimensional vectorial Sturm-Liouville eigenvalue problem {P,A,B,A,B) over 
< a; < TT, which is isospectral to the given A^-dimensional vectorial Sturm- 
Liouville eigenvalue problem (P, A, B, A, B) over < x < vr. For simplicity we 
shall assume the completeness of the eigenfunctions of the the given eigenvalue 
problem {P,A,B,A,B) subject to the given boundary conditions shown in (|I|). 

This paper is organized as follows. In section 2 we present some preliminary 
results related to vectorial Sturm-Liouville equations, and a matrix wave equation 
(see (|l8|)) which is constructed from the given eigenvalue problem (P, A, B, A, B). 
In section 3, using matrix wave equation introduced in section 2, we construct 
eigenvalue problems which are isospectral to {P,A,B,A,B). In section 4 we 
present two examples using our construction method. As shown in one of our 
examples, even for a simple case such as (P, /, 0,/, 0), where P(x) is a constant 
two by two diagonal matrix, / is the two by two identity matrix, an isospectral 
problem (Q,1, 0, /, 0) can be found where Q{x) is not simultaneously diagonal- 
izable. The isospectral problem for vectorial Sturm-Liouville equations is much 
more complicated than its scalar counterparts. 

2 Preliminary 

To study the eigenvalue problem (|^) we introduce the following matrix differential 
equation 

(4) -Y" + P{x)Y = XY, y(0) = P*, Y'{0) = -A*. 

Let Y{x;\) denote the N x N matrix- valued solution of the initial value problem 
(I). We have ( see § ), 

X 



Y{x;X)=C{x;n)+ f K{x,t)C{t; i2)dt 
Jo 



where /i^ = A , C{x;fi) = cos{fix)B* — /i ^ sm{fj,x)A* , and K{x,t) is as that 
described in [2, Lemma 2.1]. Define the following matrix- valued function 

(5) W{\) = BY'{tt-\)+AY{t:;\). 

Then A* S a{P,A,B,A^B) if and only if W{\tf) is a singular matrix. It follows 
from the variational principle that the set of the zeros of the equation 

(6) det W{X) = 



is bounded below. Denote the distint zeros of (^ in stricting ascending order as 

(7) Ao < Ai < A2 < • • • . 

Then the multiphcity ruk of A^ in the sequence of eigenvalues S(P, A, B, A, B) of 
(H) is equal to the dimemsion of the null space Null(W^(Afc)) of W{Xk). If v is a 
nonzero element in Null(M^(Afc)), then the vector-valued function 

z{x) = Y{x;Xk)v 

is an eigenfunction of (|^) corresponding to the eigenvalue A^. In the following, 
for V and w in R , the notation (v,w) is used to denote the inner product of 
two elements v and w in R^. We shall need the following result. 

Lemma 2.1 For each k > 0, in the null space NuU{W{Xk)) there are exactly ruk 
linearly independent constant vectors 9i[k), 1 < / < mt, such that the vector- 
valued functions 

Y{x;Xk)ei{k), l<l<mk, 

are mutually orthogonal, i.e., 

\Y{x;Xk)9^{k),Y{x;Xk)9j{k))dx = 0, if i ^ j. 
/o 

Proof. Let t;i, . . . , Vm,, be a basis of Null(VK(Afc)) and V^ = [vi, . . . , Vmk]- Then 

V= rVk*Y*{x;Xk)Y{x;Xk)Vkdx 
Jo 

is an rrifc x nik positive definite matrix. There exists an m^ x m^ orthogonal matrix 
U which diagonalizes V, i.e., U*VU is a diagonal matrix. Let 9i{k), 1 < Z < m^, 
denote the column vectors of VU . Then 9i{k) fulfills our requirement. □ 



According to Lemma |2.1j , we define the following vector-valued function 

(8) ^i{x;X,Xk) = Y{x-X)9i{k), 1 < / < m^. 

Then, the functions 

0/(2;; Afc, Afc) = Y{x] Xk)9i{k), I < I < ruk, 

form an orthogonal basis of the eigenspace corresponding to the eigenvalue A^ . 
From now on the eigenvalue problem (P, A, B, A, B) shall be fixed. And, as it was 
mentioned in the introduction, for simplicity, we shall assume the completeness 
of the system of eigenfunctions {(t)i{x; Xk, Xk) : I < I < mk, k = 0,1,2, .. .} subject 
to the given boundary conditions. 



Next we extend the idea used by Jodeit and Levitan in |^] to construct the 
kernel function for a related integral equation. We shall view an R -vector v 
as an A^ X 1 matrix. Choose df. € R, 1 < i < m^, k = 0, 1, . . ., which convergs 
so rapidly to zero that the matrix-valued function T , defined by the following 
uniformly convergent series, 

oo rvife 

(9) :f{x,y) = ^^4(/>j(x;Afe,Afe)(/)*(y;Afc,Afc), 

k=Oi=l 

is continuous and has continuous first and second order derivatives. Then we 
construct the following integral equation : 

(10) K{x, y) + T{x, y)+ I K{x, t)T{t, y)dt = 0, Q<y <x<-k. 

Jo 

We note that if we choose the sequence {c\ : 1 < i < m^, A; = 0,l,2,...) such 
that c|, = for all k > ko, i = 1,2,..., ruk, where ko is a fixed index, then the 
series in the right hand side of (^) is a finite series, and the equation (|^) makes 
sense. This choice of {c\.) shall be used in section 4 to construct some concrete 
examples. 

The existence of the solution of (|l^) can be easily proven by using iteration 
method. On the other hand, when we choose suitably those real numbers c\,l < 
i ^ iTT-k , k > 0, we may prove the following uniqueness theorem. 

Theorem 2.2 Suppose that the sequence {c\) is chosen so that the series in ^ 
is uniformly convergent and has continuous first and second order derivatives, 
and 

(11) i + 4\M-,Xk,Xk)\\^>o, i<i<mk, yk>o. 

Then ([J^j has a unique solution for every x , < x < vr. 



Proof. The method for proving this theorem is similar to the one used for 
treating the scalar case in [ 3, Theorem 1.1] . It suffices to show that the only 
solution for the integral equation 

(12) A(x, y) + r A(x, t)T{t, y)dt = 

JO 

is A(x,y) = 0, where A(x,?/) is the difference of two solution of (p^. Denote 
<t>i,k{x) = (t>i{x; Afc, Afc) for convenience. Owing to the assumption about the com- 
pleteness of eigenfunctions of (P, A,B,A,B), and the orthogonality (Lemma |2.1| ) 
of the eigenfunctions 4>i^k{x), we have 

(13) A-(x,,) = f;£|tM(r^A.(,.„rf,). 

r^T^. im.kW Jo im.kW 



k=Oi 



By (12), we have 

A*{x,y)+ rT*it,y)A*{x,t)dt = 0, 
Jo 

(14) A{x,y)A*{x,y)+ f A{x,y)T*{t,y)A*{x,t)dt = 0. 

Jo 

Integrating (|l4|) with respect to y- variable from to x, using (^) and (p^), we 

obtain 

(15) ^^-^[i + 4ll</..,fcl|2]a,fe(x)o*,(x) = 0, 

k=oi=o ll'^«.'=ll 

where 

(16) ^i.k{x)= r A{x,t)^i^k{t)dt. 

Jo 



Since ^i^k{x)^* i^{x) is nonnegative definite, (|T5|) and (|Tl|) imply that 

0,,fc(x)J^*fc(x) = 0, 
and hence Qi^kix) = 0, and by (|l6|). 



(17) / A{x,t)<t>i,k{t)dt = 

Jo 

for k = 0,l,2,...,i = l,2,..., rrik. Then by the completeness of eigenf unctions 
of {P,A,B,A,B), (l^ implies A{x,y) = 0. This completes the proof. □ 

Now we face the question : " Does the matrix- valued function K{x,y) deter- 
mined by the above theorem also satisfy some wave equation with which we are 
familiar as in the scalar case ? " The answer is affirmative, as shown below. 



Theorem 2.3 Assumption as Theorem 2.2. The solution K{x,y) of ( f?7| j sat- 



isfies the following partial differential equation 

(18) ^K-Qix)K = ^K-KPiy), 

where Q{x) = P{x) + 2d/dxK{x,x), and it also satisfies the following conditions: 

K{x,y) = 0, if y>x, 

(19) K{x, 0)A* + ^K\y=oB* = 0, 

dy 

(20) K{x, x) = \ r[Q{t) - P{t)\dt - ^(0, 0), 

2 Jo 

where 

oo m(j 

(21) n^,Q) = B*{Y,Y.^,e.mem)B. 



k=Oi=l 



Proof. Denote 

J(x, y) = K{x, y) + T{x, y) + T K{x, t)J^{t, y)dt. 

Jo 

Then by (p^), ^ = 0, hence Jxx = Jyy = 0. On the other hand, as 

Jxx = ^K+[P{x) + {—K{x,x) + —K{x,x))]T{x,y) 

CO m^: 
- X! X! ^k4:(piix; \k,\k)<l)*i{x] Afc, Afc) 
fc=Oi=l 



^vv = g^K + iT{x,y)+£ K{x,t)J^{t,y)dt)Piy) 

oo rrik 

+ X] 51 '^fc4'/'i(2;; Afe, Afc)0*(x; Afc, Afc) 

oo mj; 



k=Oi=l 

we have 



= J„-J„ + JP[y) 

9 „ p 52 



+ / if(x,t)[^^4(Afc0»(t;Afc,Afe)),/.*(y;Afc,Afc)]di. 

Replacing Afc0i(t; Afc, A^) by -(t)'l{t; Xk, Afc) + P{t)(t)i{t; Xk, Xk) in the last integral 
and using integration by parts twice, we obtain 

= Jxx-Jyy + JP{y) 

d"^ d"^ d 

= ^^-^K + {P{x) + 2^):F{x,y) 
ox'^ oy^ dx 

+ l^ig^Kix, t) - g^Kix, t) + K{x, t)P{t)]Ht, y)dt. 



Finally, let Q{x) = P{x) + 2d/dxK{x,x), we have 

= Jxx-Jyy + JP{y)-Q{x)J 

= + [K{x,0)—Hx,y)\x=o - -g^Kix,y)\y=oHO,y)] 
(22) + i-^K - ^K - Q{x)K + KP{y)] 

,x Q2 g 

where the function 



+ j [g^K - -^K - Q{x)K + KP{t)]T{t, y)dt, 



d d 

K{x,Q)—T{x,y%=Q - -^K{x,y)\y=oT{0,y) 



vanishes if and only if ( |19D holds. Then (|22D becomes an integral equation of the 
same type as ( p!2[). Hence we have (18). 
By (^) and (|10D, we see that 

oo rrife 
k=Oi=l 

(pO|) is a consequence of the fundamental theorem of calculus from the definition 
of Q{x) given above. □ 



Theorem 2.4 If K{x,y) is determined by Theorem 2.3 , then, for every com- 
plex X, k > and 1 < / < m^, the vector-valued function ^;(x;A,Afc) defined 
by 

px 

(23) iJi{x] X, Xk) = Mx'A, Xk) + K{x,t)(j)i{t;X,Xk)dt 

Jo 

is a solution of the vectorial differential equation 

(24) - ip" + Q{x)ip = X^j, < X < vr, 

where Q{x) = P{x) + 2d/dxK{x,x), and ipi{x; A^, A^), . . . , tpmkix', A^, A^) are lin- 
early independent. In addition, it also satisfies the following initial conditions: 

(25) V'KO;A,Afc) = B*ei{k), 

(26) i''i{0;X,Xk) = {-A*+K{0,0)B*)ei{k) 

oo rrir 

= {-A*-B\Y.Y.^A{r)eKr))BB*)e,{k), 

r=Oi=l 

or, equivalently, 

(27) B4;'i{0;X,Xk) + Ai;i{0;X,Xk) = 0, 

where 

(28) A = A-BK{0,0). 



8 



Proof. By (22), we have 



Jo 
(29) = iPix)-\I)Ux;X,Xk)+f^K{x,t)Mt-AAk)dt 

d 
+ K{x, x)(l)'i{x] A, Afc) + -rT^i^^ t)\t=x(l}i{x\ A, Afc), 

and 

/■x 

Xipii^x; X, Xk) = X(j)i{x; X, Xk) + K{x,t)X(j)i{t;X,Xk)dt 

Jo 

= X(l)i{x; X, Xk) + r K{x,t)XY{t;X)ei{k)dt 
Jo 

= A0;(x;A,Afc)+ r K{x,t)[-Y"{t;X)+P{t)Y{t;X)]e,{k)dt 
Jo 

= X(t)i{x; A, Xk) + r K{x, t)P{t)^i{t; A, Xk)dt 
Jo 

K{x,t)Y"{t;X)ei{k)dt. 



Using integration by parts twice on the last integral, we have 



AV'K2;;A,Afc) = A0/(x; A,Afc) + / K{x,t)P{t)(t)i{t;X,Xk)dt 

Jo 

(30) + {-K{x, 0)A* - ^K\y=oB*)ei{k) 

dy 

d 

K{x, x)<l)'i{x; A, Xk) + -^K\y=x(piix; A, Xk) 



dy 

rx Q2 

-J -^K{x,t)(t)i{t]X,Xk)dt. 

Then, by Theorem |2.3| , we have 

-^'{{x; A, Xk) + {XI - Q{x))i}i{x; A, Xk) 

= {-K{x,Q)A* - ^K\y=oB*)0i{k) 

rx Q2 q2 

= 0. 

Besides, for 1 < / < ruk^ 

Mo-A,Xk) = MO;\Xk) = B*ei{k), 

9 



4{0;X,Xk) = </>l(O;A,Afc) + K(O,O)0KO;A,Afe) 
= {-A* + K{0,0)B*)ei{k) 

oo rrir 

= [-A* -B*{Y,Y.^A{r)ei{r))BB*)em- 

r=Oi=l 

If we denote A = A- BK{0, 0), then BA* = AB*, and, by using BA* = AB*, 
we have 

BV'KO;A,Afe)+iV'z(0;A,Afc) 

= B{-A* + K{0, 0)B*)ei{k) + {A- BK{0, 0))B*ei{k) 
= 0. 

The hnear independence of those vector- valued functions ipi{x; Xk,Xk),l < I < m^ 
can be proven by using (p3[), Gronwall's lemma and the linear independence of 
those functions (/)z(x; Afc, Afc), 1 < / < mfc. □ 

At the end of this section, we state a theorem which indicates the possi- 
ble candidates for eigenfunctions of the isospectral problem {Q,A,B,A,B) of 
{P,A,B,A,B) which shall be described in next section. Furthermore, we may 
also use this theorem to construct that an isospectral Dirichlet's problem of a 
given Dirichlet's problem. 

Theorem 2.5 Suppose Xk G a{P,A,B,A,B), k>0. Then 

(31) ipi{x;Xk,Xk) = (j)i{x;Xk,Xk) 

oo nir 1.x 

- ^ ^ 4Vj(a:; Xr, Xr) / (jfiit; Xr, Xr)(l)l{t] Xk, Xk)dt 

for all 1 < I < m-fc. 

Proof. The proof is similar to the one of [3, Theorem 1.4]. Denote 4'i,k{x) = 
4>i{x; Xk, Xk). By (^), and the integral equation (10), we have 



K{x,y) = -T{x,y)- f K{x,t)T{t,y)dt 

Jo 

oo rUr px 

^^4[0.,.(x)+ / Kix,t)(i)iAt)dm,riy)- 

irn"; JO 



r=0 i=l 

Then, by (|23|), we have 



(32) K{x,t) = -^g4V'*(x;A„A,)0f(t;A„A,). 



r=Oi=l 

Apply m, m implies m. D 
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3 Isospectral problem 

Those theorems in previous section enable us to construct an isospectral problem 
from a given eigenvalue problem (P,A,B,A,B) and a sequence of real numbers 
cl.,1 < i < ruk, k > 0, where the sequence (c|,) satisfies the assumption of Theo- 



rem |2.2|. As Theorem 2.4 states, for any A; > 0, and for each I, 1 < I < rrik 



the vector- valued function il^i^x; A^, A^) satisfies the boundary condition 

5V''(0) + iv(o) =0, 

where A is given by (p3)- Hence, the final step for constructing isospectral problem 
is to determine the form of boundary condition to be satisfied at x = vr. For this 
purpose we use formula (pl|). By (^TI), we have 



l + cl\\(j)i{x;Xk,Xk)\r 
Differentiating pT| ) with respect to x and evaluating it at vr, we have 

ip'l{Tr;Xk,Xk) = (l)'i{TT;Xk,Xk) -ciip'i{7r;Xk,Xk)\\(l)i{x;Xk,Xk)\\'^ 

OO rUr 

- [JUI 4^i(7r; Xr, Ar.)(/>r(7r; Ar-, Ar-)](/>/(vr; Afc, Afc), 

r=Oi=l 

and, hence we have 

(1 + ci\\cPi{x; Xk, Afc)||2)V'Kvr; A^, A^) 
= 0z(vr;A.,A,)-[^J^^^^--^-— — ^]0,(7r;A,,A,) 



r=0 i=l 



Acting on the above identity by B and using the condition S(/);(7r; Afc, A^) + 
^i?i>i(7r; Afc, Afc) = 0,we have 

^(l + ci||</.Kx;Afc,Afc)||2)V'K7r;Afe,Afe) 

R^'^ \ \ ^ gr'^^ c;'/'i(7r;A^,A^)(/)*(7r;A^,A^) 
= ^^.(vr, A„ A.) - ^Eg l + 4||^.(.;A.,A.)|P l^'^"' "'=' "'^^ 

^ /I I Kr-^v^ 4'A»(^;Ar,A^)(/)*(7r;A,.,Ar) 
Then, by (|33|), we have 

nil I \ \\ I A \ KfX^Xr^ ^r't'ii^^'j K, Xr)(l)*iT^',Xr,Xr)-.. r \ \ \ 

S^,(vr, A„ A,) = -{A + e[g g i + cUmx;K,XrW "^'^ ' '^' '^' 

11 



and hence, 

SV'K^; ^k, Afc) + AiPi{tt; Afc, Afc) = 0, , 
where 



(34) A = A + B[Y^Y.-^ 



r=oi=i ^ + 4\\(l)iix;K,KW 

In fact, by (|32|), (^^ can be simphfied as 

(35) A = A-13K{tt,7t). 

Furthermore, if we can prove that, for any R -valued function / satisfying 
Bf'iO) + i/(0) = , Bf'iTT) + Af{7T) = 0, and 

{f{x),ipi{x;Xk,Xk))dx = 0, k>0, 1 = 1,... ,mk, 



we have / = 0, then the set {tpi{x; A^, A^) : k > 0,1 < I < mi.} is complete, and. 



by Theorem g]4|, we have T.{P, A, B, A, B) = T,{Q, A, B, A, B). For our purpose, 
let 

nf){^)= r K{x,t)f{t)dt. 

Jo 
Then 



T*{g){t)= j^ K*{x,t)g{x)dx. 
Writing ipiix; Xk, Xk) = (I + T)(j)i{x; Xk, Xk) , 1 < I < nik, we have 

0= {f{x),ipi{x;Xk,Xk))dx = / {f{x),{I + T)(t)i{x;Xk,Xk))dx 
Jo Jo 

{{I + T*)f{x),Mx;Xk,Xk))dx. 



Now set g = {I+T*)f. Ifwe show that ^^'(0) + ^5(0) = 0, and Bg'{TT) + Ag{TT) = 
0, then by the completeness of the set {(l)i{x; Xk, Xk) : k > 0,1 < I < mk}, we 
have 5 = and hence / = 0. We only check the identity Bg'{0) + Ag{0) = 0, 
the other part can be proved by similar argument. Suppose Bf'iO) + AfiO) = 0. 
Then, using (p^) and (p8|), we have 

Bg'(0)+Ag{0) = B[f'{0) - K*{0,0)f{0) 



+ K;{x,0)f{x)dx]+A[f{0)+ K*{x,0)f{x)dx] 
Jo Jo 

= Bf{0)-BK*{0,0)f{0) 

+ ^/(0)+ r [BKl{x,0) + AK*{x,0)]f{x)dx 
Jo 

= Bf'iO) + [A- BK{0, 0)]/(0) = Bf'iO) + Af{0) = 0. 
As a conclusion of the previous arguments, we have the following theorem. 

12 



Theorem 3.1 Let mt denote the multiplicity of Xk in a{P,A,B,A,B). Suppose 
{cfc, 1 < i < rnk,k > 0} is a sequence, satisfying the condition | fjj| j and making 
J-{x,y) in ^) a C"^ -function, Q{x) is as that defined in Theorem 2.3, A and A 
are as those defined in §^) and B^). Then T.{P, A, B , A, B) = ^{Q,A,B,A,B). 



As a final remark, we note that if A = I , B = 0, A = I , and ;S = in (|T|), then 
the matrices A and A in Theorem |3.lJ are equal to I, the identity matrix. Hence, 
for a given Dirichlet problem (P, /, 0, /, 0), the isospectral problem constructed in 
Theorem 13. 11 is also a Dirichlet problem. 



4 Examples 

In this section, we use our theory to construct some examples which have some 
significant meaning the scalar case can not tell. 

Suppose Ao is an eigenvalue of (jl]) with multiplicity rrio. Let (/>o(x) = col{4>i{x), 
(/)2(x), • • • , (J)n{x)) be an eigenfunction corresponding to Aq. Take 

T{x,y) = c4>o{x)4>l{y) 

( (t)i{x)(t)i{y) 0i(x)02(y) ••• (t)i{x)(j)N{y) \ 
(l>2{x)cl)i{y) 4'2{x)4'2{y) ••• (/>2(a;)</'7v(y) 

\ cl)N{x)cl)i{y) (/>Af(a;)(/>2(y) ••• cl)N{x)(l)N{y) ) 

Plugging it into (p^), and letting kij{x,y) denote the {i,j) entry of K{x,y), 
we have 

(36) kij{x,y) + c(pi{x)(pj{y) + (c / (^kir{x,t)4>rit))dt)(j)j{y) = 0, 



r=l 



for i = 1, . . . ,N and j = 1, . . . ,N. 
We shall show that 

(37) k.,ix,y)- '^^^^^^^^^^ 



l + cJ,^\Mt)\'df 
For i fixed, consider the equations 

N 



(38) kij{x,y) + ccl)i{x) 



hj{y) + {c [ (V k„{x,t)Mt))dt)<pj{y) = 0, 



^ 1^ j ^ N. Multiplying the jth equation by (pj{y), integrating it from to x with 
respect to y, and denoting Pij{x) = f^ kij{x,t)(f)j{t)dt and aj{x) = Jq (l)'j{t)dt, 
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^ ^ j 1^ N, we have the following lineaar system of equations with unknowns 



N 



(39) Pij{x) + c(j)i{x)aj{x) + caj{x)(^pir{x)) = 0, 1 < j < iV. 

r=l 

Solving this system, we obtain 

, X / s C(hj(x)aA(x) 

Note that, by (|8|), 

kij{x,y) = -c(j)i{x)(j)j{y) - {c (V kir{x,t)(j)r{t))dt)(l)j{y). 
Hence if we plug ( |40| ) into (|3^) , we obtain (|3^) . It also follows from (^) , that 



Hence, according to Theorem |3.1| , by setting 

Q{x) = P{x)-2^[- '^°(-)'^°(^) 



(41) A = A-cBM^WoiO), 

c</>o(vr)(/'o(vr) 



A = A 



l + c! 



we have S(P, A, B, A, B) = S(Q, A, B, A, B). 

As an example to the above construction, we construct the following eigenvalue 
problem which has an eigenvalue of multiplicity 2. 

Let I be the 2x2 identity matrix. Take 

P{x) = i~^ Q j , A = A = h B = B = 0. 

Then one can verify that for the eigenvalue problem (P, /, 0,/, 0), 1 is an eigen- 
value of multiplicity 2, and the other eigenvalues are all simple. The eigenspace 
corresponding to 1 is the vector space spanned by the two vector-valued func- 
tions (sin(2x),0)* and (0, sin(x))*. Choosing (sin(2x),sin(x))* as the eigenfunc- 
tion which plays the role of (j)o{x) in the above construction, c = 1, and using 
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(p]), we have 

QM = (-^ M - -[ - 

Csin^(2x) sin(2x) sin(x) 
sin(x) sm(2x) sin^(x) 

and A = A = I. Note that the matrix potential function Q{x) is not simulta- 
neously diagonalizable since, as checked by computation, the matrix-valued func- 
tions Q{x) and Q'{x) do not commute. On the other hand, if we take (0,sin(x))* 
instead, then we find that Q{x) is a diagonal matrix- valued function and is of the 
following form 

/ -3 \ 

Q{x) = 

n d ( -2sin2(x) X 

\ ^ a!^4+j;;sin2{t)dt'' / 

There are lots of interesting phenomena can be observed from our construction, 
which shall be observed later. 
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